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Abstract. Using the Inverse Scattering Method with a nonvanishing boundary condition, we 
r- H obtain the square k 2 of a focusing modified Korteweg-de Vries (mKdV) breather solution with 

—i non zero vacuum parameter b 2 . We are able to factorize and simplify it in order to get explicitly 

the associated mKdV breather solution k with non zero vacuum parameter b. Moreover, taking 
^SJ the limiting case of zero frequency, we obtain a generalization of the Double Pole solution 

I i introduced by M.Wadati et al. 

o 
O 

1. Introduction 

The focusing modified Korteweg-de Vries equation (mKdV for short) 

dk dk , ndk m m , 

^ m + d^ + 6k d- s =° As > t)eRxR - (L1) 

^3 appears to be relevant in a number of different physical systems (e.g. phonons in anharmonic 

lattices, models of traffic congestion, curve motion and fluid mechanics between others). Indeed, 
it can also be considered as a canonical equation, as KdV, sine-Gordon and non linear Schrodinger 
equations. 

Breather solutions of mKdV ( TT7TT ) in the line were found by M.Wadati [T3] (see also G.Lamb 
^ I [TU]). They were rediscovered by C.Kenig, G. Ponce and L.Vega in the proof of the discontinuity 

£f~) of the flowmaps associated to mKdV equation in the Sobolev spaces H a , constituted by functions 

Q\ with a derivatives in L 2 (M) (see [9]). Those breather solutions are defined in all the real line, vanish 

exponentially at infinity and, in a qualitative point of view, describe traveling wave packets. 

Indeed, they are determined by four real parameters, two of them given by the amplitude ((3) 
of the envelope and the frequency (a) of the carried wave and the other two given by time and 
spatial translations (represented by uq , «o below) . It is well known that such breather solutions 
can be written as follows 

, / s 9 , ff(u) + ig(v)\ „<9 fg(v)\ , 

k(s,t ) = -i 7r log ^ ^4 =2— arctan , (1.2) 

ds \f(u)-ig(v)J ds \f{u)J 



X 
S3 



where u = 2f3s + jt + u , v = 2as + 6t + v , 7 = 8/3(-/3 2 + 3a 2 ), 6 = 8a(a 2 - 3/3 2 ), ancQ 



g(v) = — sin I v — arctan I — I I , (1.3) 



a 

f(u) = cosh(«). (1.4) 

This paper is devoted to the use of the Inverse Scattering Method (ISM for short) under a 
nonvanishing boundary condition(NVBC shortly) devised by T.Kawata and H.Inoue j8] to obtain 



breather solutions of mKdV (1.1 ) which at infinity behave as a non trivial constant parameter b. 
The interest of this kind of solutions of the mKdV is related with the problem of the evolution of 
closed planar curves under the mKdV flow in the following way. By the work of R.Goldstein and 



D.Pctrich [B] the equation ( 1.1 ) is considered as the evolution equation of the curvature of a curve. 



Date: September, 2011. 

2000 Mathematics Subject Classification. Primary 37K15, 35Q53; Secondary 35Q51, 37K10. 
Key words and phrases. mKdV equation, Breather, inverse scattering, nonvanishing condition. 
^The arctan phase in the argument of g can be dropped with a suitable translation in time and space, but it 
will not be done for comparison purposes. 

1 



2 



Focusing mKdV breather solutions with nonvanishing boundary conditions by the ISM 



By this work, closed curves can be considered as those whose associated curvature satisfies that 
its mean is non-zero. Moreover, those curvatures can be obtained from a solution of the focusing 
mKdV constructed by the (nonlinear) superposition of a constant (e.g. constant b) plus a travel- 
ing wave. In fact, it is posible to find some special breather solutions associated to simple closed 
curves that, when evolving under the mKdV flow, they create and annihilate self- intersections (see 
[5] for more information). 

Although in the literature some kind of breather solutions of the Gardner equation(also known 
as the extended KdV equation) have been obtained before [JJ H2] , they do not contain details 
on derivations of these breather solutions from the ISM scheme. Even despite the close relation 
between Gardner and mKdV equations (solutions of the mKdV equation with NVBC are also 
solutions of the Gardner equation with zero boundary condition), breathers obtained in these 
references [7J [T2] are far to be easily compared with the breather of M.Wadati [13] they gener- 



alize. These two aspects are attained too in the present paper(for further details see section 2.2 



and (2.47)). Finally note that the mKdV breather solutions mentioned above can be obtained 
alternatively using the Hirota method with a suitable selection of the wavenumbers (see the work 
of K.Chow and D.Lai 0). 

2. Breather solutions of the focusing mKdV with nonvanishing boundary 

conditions 

In this section we obtain breathers of the focusing mKdV ( |1.1[ ) with nonvanishing boundary 
conditions by using the ISM for potentials that are not trivial at infinity, introduced by T.Kawata 
and H.Inoue in [8]. We also recall the work of T.Au-Yeung et aL[4], in which the same approach 
was used to obtain one and two soliton solutions with non trivial values at infinity. First, we 
summarize some basic results from [5] and [1], necessary for our research. 

2.1. Basic results of T.Kawata and H.Inoue for the mKdV. T.Kawata and H.Inoue con- 
sidered in pQ a generalized AKNS eigenvalue problem for nonvanishing potentials, which consist 
in the following spatial and time evolution equation^ 

v s = D(X)v,D(X;s,t)=(-£ ) (2.1) 
v t = F ( A)v,F(A;M)=(^ , (2.2) 

where 

A(X) = -4iX 3 - 2i\qr + rq s - qr s , 

B(X) = 4A 2 q + 2i\q s + 2q 2 r - q ss , (2.3) 

C(A) = 4A 2 r - 2iXr + 2qr 2 - r ss . 
Matrices -D(A) and F(X) satisfy the well known integrable condition associated to equations 



( pTj ) and Q (i.e <9 t ([2T]) = d s $Lty): 

D t — F s + DF — FD = 0. (2.4) 

They seek a real potential solution q(s,t) under the following boundary condition: 

q(s,t) ->6ass-> ±oo; b 2 = -A 2 ,, (2.5) 

requiring that q(s, t) is sufficiently smooth and all the s derivatives of q tend to zero as s — > 
±oo. For this purpose, they consider potentials q(s,t) and r(s,t) with the following nonvanishing 
conditions: 

4n what follows, subindiccs of the form f s means d a f = 2^. 
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q(s,t)(oT r(s,t)) — > q ± (ov r ) as s — > ±00, 

(2.6) 

q+r+ = q r = A§, 

where q^,^ and Aq are constants. Then, the spatial evolution matrix D(X;s,t) can be written 
as follows: 

D(\;s,t) = D ± {\) + AD ± {s,t), (2.7a) 



D ± (X) = lim D(X;s,t) 

S— >±oo 



-iX q ± 

-± iX 



(2.7b) 



A£» ± (s,t) = D(X;s,t) - D ± (X) = 





r(s, t) — r 



q(s,t) - <? d 




and the characteristic roots of D ± (X) are ±i£> with £ = y/X 2 — Aq. Now, they define 

T ± (A, C) = lim s ^±oo T(A, C; s), 

riA.^si-'-^W A ^ 



(2.7c) 



(2- 



where 



A-C in(«),/' 



and (71, ri are suitable smooth and satisfy that qi(s)ri(s) — Aq for all s € R. The matrices D ± (A) 
can be diagonalized by T ± (A, C) as 



^ ± (A) = - l cr ± (A,c)tT3[r ± (A,c)]- 1 . 

Using (2.9|, they can define Jost matrices < i> ± as the solutions of (|2.1| under conditions 



$ ± (A,C;s,i) -> T^O-ACs) as s -> ±00, 



(2.9) 
(2.10) 



where 



J((s) 







Then, a scattering matrix S is defined by 

$-(A,C;M) = $ + (A,C;M)S(A,C;*), 



and using relations (2. 2), (2. 9), (2.11) it is easy to derive the following condition: 



St + sw~ - w+s = 0, s\ t=0 = s , 

where So is given by the direct scattering, and 



(<t 3 • ^ a„A™ 1 + i{a Q /C) ( ± e -2iCs 

n—1 ^ 



iX 



The solution of (2.13) is easily obtained as follows 



S(XX;t) = e w+t S (X,C)e 



-w~t 



(2.11) 
(2.12) 

(2.13) 

(2.14) 



At this point, we apply this framework to the mKdV equation (1.1), by choosing q(s,t) — k(s,t) 
and r(s,t) = —q(s,t). Then, we calculate from (2.13) the temporal evolution of the elements of 
the scattering matrix (SV,). For that purpose (see [5J p. 1723]), we select N = 3 and do = a 2 = 0, 
ai = 2ib 2 , 03 = — Ai in (2.13), so that: 



Sii(A,C;f) = S n (A,C;0), S ia (A,C;*) - S 12 (A,C;0) e - 4 ^ 2A2 - h2 >', 
S 21 (A, C; t) = S 21 (X, C; Q)e u ^~ b ^\ S 22 (X, (; t) = S 22 (A, C; 0). 
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Now, assume that zeros of the S l n(A, £; 0) matrix element in the region 3(C) > are (Xj, Q), j 
1 , 2 where 



= ^A2 + 6 2 , j = 1,2, forj/> a , (2.15) 

and define the functions K ± (s,y) as those satisfying the following three equations (for the shake 
simplicity, we drop the time dependence in K and H c defined below): 



2§± (s, y) + a 3 2*± (s, y)a 3 + a 3 K± (s, y)a 3 [D± (A) + iA<7 3 ] - ( 



fc(s,f) 
fe(s,t) 



(73/^ (s, s)a 3 - AT* (s, s) + AD* (A; s, t) = 0, 
if (s, y) — > as y — > ±oo, 

'1 



where D ± (X) and AD ± (A;s,t) are defined by (2.7 1 and 03 
Levitan equation reads as follows, 

K + (s, y) ■ ( ? ) + E?=i ( e- ) ^ (s+1/) - ff c (« + ») ( ? ) 

r +oo T ^+, _.i\\-^2 / c j \ Afu+i/lj.j , r+°° 



-1 



*T±( S ,y)=0, 

(2.16a) 
(2.16b) 

(2.16c) 

Then, the Gelfand- 



where 



K+ (s, y>) E - =1 ( ) ^ W + /+ 00 K + (a, y')H c (s + y>) ( ? ) <V = 0, 



i( Aj -g)m(A i )e^^-^ t ! 



m(A) 



^ TO (A,-)e 4iC ^ 2A ?- 62 )*, 

S 2 i(A,C;0) 
C^5 n (A,C;0)" 



(2.17) 

(2.18) 
(2.19) 
(2.20) 



We rewrite (2.171 assuming: 



a. ) Continuous component 

ff c (y;0) = 0. (2.21) 

b. ) Zeros of 5n(A, C) m the region 3(C) > are (Xj, Q) and (— Xj, Q), j = 1, 2, where 



Ci = \/Af + & 2 , Ai = a + i/3, A 2 = -A* = -a + ^. 



(2.22) 



Taking into account that m(— A*) = m(A)* (see p. 1728]), (2.17) becomes 



4i< (2A^-b^)t 



4iC,(2Af-i>" i )t 



-iC j m(Xj)e*"'J l "i 



e iCj(s+3/) 



iCi(y+y')d y ' = 0. 



(2.23) 



Choosing a representation of I J as follows 



(KUs lV ) 
\Kt 2 {s,y) 



2 

E 



then 



+ , . (K+(s, y) K+(s, y)\ ( K+(s, y) K+(s, y) 
[,V) \K+( S ,y) K+( S ,y)J \-K+( S ,y) K+(s,y) 

Now, substituting equations (12.241) and (12. 251) in (12. 231), we obtain the system 



(2.24) 



(2.25) 
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Kj(s) + ELl ^"( S ) a i ^i(Cj+Cn)' ~ En=l #"( S ) a . 



■3 *(C 3 +C„) 



+ a,-e^ s = 0, 



where 

which can be rewritten in a matricial form as 



where K = 



K 2 



E B\ I K 



K = 



K-2 



A = 



-a\e 
-a 2 e 



<2S 



A = 



ids 



a 2 e 



(2.26) 

(2.27) 
(2.28) 



E — (Ej n ) — Sj, 



Defining the determinant of the coefficient matrix by 

E B 



A 



—B E 



and by using (2.24), we obtain that the solution of the system (2.26) is 

= -~0°8 A). 



(2.29) 
(2.30) 

(2.31) 
(2.32) 



From the components of (2.16) we get, k 2 {s,t) = b 2 — 2^^22(5, s). Inserting (2.32) into this 
relation, we finally obtain the expression 

,2 , d 2 (logA) 



k 2 (s,t) 



ds 2 



(2.33) 



2.2. The breather solution of mKdV with nonvanishing boundary value. Now, our aim 
is to obtain an explicit expression of the focusing mKdV solution k from (2.33). For that purpose, 
we rewrite the determinant (2.31 ) in the four dimensional case (as it corresponds to the breather 
case) as a product of two simpler 2x2 determinants. Before that, we remark the following facts 
about the roots (Ci,C2) and the temporal dependence of the coefficient Sn- First, recall that in 
the breather case with nonvanishing boundary conditions the roots of Sii(A, £) are given by 



Cj = ^/\j + b 2 , j = l,2, A 1= a + z/3, A 2 = -A* = -a + i/3. 
We express Q, i = 1,2 as complex numbers as 

Ci = a(a,/3,b) +0(a,/3,b), 

C2 = -a(a,P,b) + ij3(a,0,b), 
where we only consider the region 3(C) > so that 



a(a,0,b) = i/{a 2 - f3 2 + b 2 ) 2 + Aa 2 f3 2 cos Q 
j3(a, (3, b) = ^(a 2 - (3 2 + b 2 ) 2 + Aa 2 f3 2 sin Q 



2a/3 

2 arctan( a2 _ /32 + fo2 



(2.34) 
(2.35) 

(2.36) 
(2.37) 



which satisfy 



a(a, (3, 0) = \/ a 2 + ft 2 cos ( arctan( — ) ] = 5R(a + i(3) = a, 

V a J 



Focusing mKdV breather solutions with nonvanishing boundary conditions by the ISM 

P(a, P, 0) = V" 2 + P 2 sin ( arctan (~)) = 3 ( a + ^) = @- 



Now, by using (2.35), the exponents of the temporal dependence given by (2.271 are rewritten as 

4id(2A 2 -b 2 ) = -%a,p,b)+iS(a,P,b), 



4< 2 (2A 2 - b 2 ) = -7(a, p, b) - i~8{a, P, 6), 
and, again, we only consider the region 3(£) > so that 

7(a, P, b) = 8/3(3a 2 - p 2 ) - 12b 2 p, 
8{a, P, b) = 8a(a 2 - 3/3 2 ) - 12fe 2 a, 

which satisfy 

7(a,/3,0)=8/3(3a 2 -/? 2 )= 7 , 
8(a,P,0) = 8a(a 2 - 3/3 2 ) = (5. 



(2.38) 

(2.39) 
(2.40) 

(2.41) 
(2.42) 



For the sake of simplicity and without loss of generality, we rename de,P,^,5 as a, /3, 7, <5, re- 
spectively. Now, we are able to rewrite the determinant A as a product of two simpler 2x2 
determinants. 



A 



= det ((I - iM - ibN) ■ (I + iM-ibN)) 
= det(I - iM - ibN) det(I + iM - ibN) 

= {det(I - iM(s, t) - ibN(s, t))} 2 + 3 {det(I - iM(s, t) - ibN{s, t))} 2 , 



(2.43) 



where 



7=1 



1 




-2/3s + ( — ,+iS)t - 



2/3 



2,2-1,, J, M(s,t) - I i{ _ a+ifj)m * e -2f»-^+iS)t - . j j i . ~ I 



N(s,t) 



2{a+i/3) 



2:1 



2/3 

-2/33 + (—, + iS)t - 

2(-a+i/3) , 



(2.44) 



Then, substituting (2.43) in (2.33), and resorting to the identity 



„ d , z . 9, w — iz 

2 — arctan - = z— log — — ), 
as w us w + iz 



(2.45) 



we get 



'-(s, t) =b 2 + £s log ($1 {det{I ~ iM(s, t) - ibN(s, t))} 2 + 3 {det(I - iM(s, t) - ibN(s, f))} 2 ) 



6 + if log 



9t{det(I-iM(s,t)-ibN(s,t))}-iS{det(I-iM(s,t)-ibN(s,t))} 
K{det(/-iM(s,t)-i6iV'(s,t))}+iS{det(/-iM(a,t)-i6JV(»,t))} 



& + 2 arctan 



Sf{rfet(J-iM(s,t)-i6Ar(s,t))} 
5R{det(/-iM(s,t)-i&Af(s,t))} 



(2.46) 

which gives directly the matricial expression for the breather solution of the focusing mKdV with 
nonvanishing boundary value: 



k(s,t) 



b + 2— arctan 

OS 



3 {det(I - iM{s, t) - ibN(s, t))} 
SR {det(I - iM(s, t) - ibN(s, t))} 



(2.47) 
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In fact, it is possible to calculate from (2.47) the explicit expression for k. For that, we first 
calculate the determinant in (2.47) (we write m = |m|e*^): 

f(s, t) + ig{s, t) = det(I - iM(s, t) - ibN(s, t)) 



1 _i_ 1 1 b 2 \ a 2 \m\ 2 -A0a-2it , m p 2i(a+i0)s+(-y+i6)t i jm* 2i(-a+i0)s-h+: 

1 ' \ L a 2 +0 2 > 40 2 e "r" 1 2 e ' 2 B 



<S)t 



(2.48) 



-4,8s— 27* i b\m\ -20s— -yt/ 

'20{a 2 +0 2 ))^ + a 2 +0 2 



>^'/";2'r^2^ e~^°~-'- + -ttkie-"*"- ' 4 (asin(2as + St + <f>) - /3cos(2as + St + </>)), 



so that 



9[ ' j ~2 + 2 6 + 2/?(a 2 + /? 2 )) 



-40s-2-yt 



f(s,t) =1 + (1 



* 2 |m' 2 
a 2 + /3 2 ' ^/J 2 



a 2 H 2 4 /? s -2 7 t 



+ 5 H e -2/3s-7t( a sin / 2aa , + ^ + ^\ _ ^ cos ( 2as + ft + 0)). 

or + p z 



By defining e 2 ^" = (1 — aHffl ) "^^ 1 ^he expression of <? and / simplifies to 

=2e- 2 ^-^-^5( S ,t), 
f(s,t) =2e- 2 ^ s -^f(s,t), 



(2.49a) 



(2.49b) 



(2.50a) 
(2.50b) 



where 



a 2 + p 2 



+ P 2 -b 2 



cos(2as + St 



bp 



a 2 + P 2 - b 2 
f(s, t) = cosh(2/3s + jt + V>) + 



(cosh (2/3s + 7* + V>) - sinh (2/3s + 7* + ip)) 
b(3 



= sin(2as + St + 4> — arctan( — )). 

a^Ja 2 + P 2 - b 2 a 



(2.51a) 
(2.51b) 



Hence, the explicit expression for the breather solution of the focusing mKdV with nonvanishing 
boundary value b (or b-breather) is: 



k(s, t) = b + 2— arctan 
Ox 



%[f(s,t)+ig(s,t)] 



= o + 2— arctan 
Ox 



<j(s,t) 



f(s,t) 



(2.52) 



_$t[f(s,t)+ig(s,t)]_ 

with f,g given by (2.51). In the formal limit b — > 0, ( |2.52 1 is reduced to the well known breather 
solution (1.2| of the focusing mKdV equation (see [15] . up to translations in time and space). 
Even more, it is possible to obtain the generalization of the Double Pole solution presented by 
M.Wadati et al. |11| with a nonvanishing boundary value at infinity. For that, do the translation 
v = v + arctan(^), where v = 2as + St + <p is the argument of the oscillatory functions in (2.51 ), 
and calculate the formal limit a — > 0. Such generalization is given by the explicit formula 



where 



G(s,t) 



F(s, t) = cosh(u) 



,/ \ , n d ( G(s,t) 

k(s, t) = b + 2— arctan _ ' 
9s \F(s,t) 



(1 - 2P(s - 6(2/3 2 + b 2 )t)) bp (cosh(u) - sinh(u)) 



(2.53) 



y/p 2 - b 2 
2bp(s - 6(2/3 2 + b 2 )t)) 
^P 2 - b 2 ~' 



P 2 - b 2 

u= P(2s-4(2P 2 +3b 2 )t). 



Focusing mKdV breather solutions with nonvanishing boundary conditions by the ISM 



Taking into account the point- wise convergence of (2.531 to b when time goes to^Ji — > +00, it is 
possible to guess its asymptotic form at the mentioned limit, which is 



k(s,t) 



2f3 2 



+ b 2 



scch(w — S + ) 



2/3 2 



Vf3 2 + b 2 



sech(w — S-), 



with 



5 ± =p\og 12/3( 



r ) ± s(4 j 8 a + 26 a )t 



(2.54) 



(2.55) 



K p 2 + b 2 ' 

The phase S± determines how evolves the distance between the soliton and the antisoliton of 



(2.53) 





0. Right: Double Pole solution (2.53) with p = l, b = 0.3 at t= 5. 



Figure 1. Left: Breather solution (2.52) with a = 7,0 = 1, b = 0.3 at t 



3. Summary and remarks 
In this paper we have obtained the breather solution of the focusing mKdV equation with 



nonvanishing boundary conditions at infinity (2.52) by using the inverse scattering method for 
potentials that are not trivial at infinity as it was devised by T.Kawata and H.Inoue in [8]. As 
far as we know, it has not been reported before a systematic work on the obtention of this kind 
of breather solutions of mKdV under the ISM. These solutions play an important role in the 
construction of closed curves with localized perturbations, which evolve under the mKdV flow 
of curves (see [2]). We have also generalized the Double Pole solution found by M.Wadati et al. 

to the case when it takes non trivial values b at infinity. We have shown that even in this 
generalization, the distance between humps grows proportionally to log(f), as the formula (2.55) 
shows. Moreover, the associated closed curve to this (Double Pole curvature) solution is a closed 
curve with two loops. These two loops enclose asymptotically the same area, they point in- and 
outward respectively the closed curve, travel in the same direction and the distance between them 
grows slowly (proportionally to log(f) as t goes to 00) (see [2]). We think that the asymptotic 
property (2.55) could be useful to check the accuracy of numerical methods (e.g. difference and 
pseudo spectral methods) for t big enough, as it was shown in [3] when b = 0. 

Finally we would like to remark that it is possible to obtain new solutions of the defocusing 
mKdV equation 

dk dk „, 9 <9fc 



dt ds 3 



6k< 



ds 



= 0, (M) G 



(3.1) 



from the focusing mKdV breather solutions in (2.52) with a special choice of their free param 
eters b,/3,a. With this goal in mind, we make the following changes of parameters in (2.52) 



'The case t — ¥ —00 is equivalent. 
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b becomes ib, f3 becomes — /3 and a becomes ia. Then, we get a new purely complex solution 



ik(s,t) of the focusing mKdV equation (1.1). Hence k{s,t) is a real and regular two soliton so- 



lution of the defocusing mKdV equation (3.1) with nonvanishing boundary value b at infinity. 
With other selections of the parameters, b, ft, a and with the same procedure, we obtain different 
complex and regular or singular (depending on the selected parameters) solutions of the equation 



(3.1) with nonvanishing boundary value b at infinity. 
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